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The first-principle theory of electron dephasing by disorder-induced two state fluctuators is devel- 
oped. There exist two mechanisms of dephasing. First, dephasing occurs due to direct transitions 
between the defect levels caused by inelastic electron-defect scattering. The second mechanism is 
due to violation of the time reversal symmetry caused by time-dependent fluctuations of the scat- 
tering potential. These fluctuations originate from an interaction between the dynamic defects and 
conduction electrons forming a thermal bath. The first contribution to the dephasing rate saturates 
as temperature decreases. The second contribution does not saturate, although its temperature de- 
pendence is rather weak, oc T 1 / 3 . The quantitative estimates based on the experimental data show 
that these mechanisms considered can explain the weak temperature dependence of the dephasing 
rate in some temperature interval. However, below some temperature dependent on the model of 
dynamic defects the dephasing rate tends rapidly to zero. The relation to earlier studies of the 
dephasing caused by the dynamical defects is discussed. 

PACS numbers: 73.20.Dx 



I. INTRODUCTION 

The problem of dephasing of electron states in low- 
dimensional structures is in focus of interests of many 
research groups. This is due to novel experiments on the 
Aharonoii-Bohm effect in specially designed mesoscopic 
circuitstra and on weak localization magnetoresistance in 
low-dimensional samplesJ3^|-jwell as to new theoretical 
discussions of dephasing BEoQoB In particular, dephas- 
ing due to defects with internal degrees of freedom as a 
source of dephasing were recently addressed. 00 Accord- 
ing to the model discussed in Ref. [?], a temperature in- 
terval can exist in which the dephasing rate is almost 
temperature-independent . 

In this work we revisit the dephasing due to dynamic 
defects which interact with electrons and tunnel between 
their two states due to interaction with some thermal 
bath. Examples pf-such defects are disorder- induced two- 
state fluctuatorgI90 present in any disordered material, 
impurities with a non-compensated spin, etc. These de- 
fects produce a random time-dependent field and in this 
way they violate the time-reversal symmetry of the prob- 
lem. According to a conventional opinion, this property 
is sufficient to produce dephasing. However, this is true 
only under the condition that a typical defect relaxation 
time is shorter that the time during which the electron 
interference pattern is formed. Indeed, if the defects do 
not change their state during the pattern formation they 
act as static ones and can gontribute to the interference 
only in a constructive way.0 

The purpose of this paper is to develop a systematic 
theory of weak localization with dephasing due to dy- 
namic defects interacting with electrons which results in 
a smooth temperature dependence at relatively low tem- 
peratures T. The dynamic defects are specified as two- 
level tunneling states (TLS) that exist in any crystalline 
metal. 



The main message of this paper is the following. There 
exist two mechanisms of electron dephasing due to dy- 
namic defects. The first one is due to direct inelastic 
transitions between the levels of the TLS leading to the 
possibility of determining the actual path of the elec- 
tron, and consequently to loss of interference. The sec- 
ond one is due to relaxation dynamics of dynamic defects 
which fluctuate due to interaction with the thermal bath. 
Time dependence of the electron scattering crossection 
due to the defects' fluctuations lead to violation of the 
time-reversal symmetry and, as a consequence, to deco- 
herence. To our knowledge, the theory relevant to the 
second mechanism has not been developed. However, 
there exists a temperature interval where this relaxation 
mechanism is dominating. 

The paper is organized as follows. Below we will give 
physical considerations to describe dephasing by dynamic 
defects which will be then confirmed by a diagrammatic 
approach, see Sec. |l[ In this section the mode l for 



electron-TLS interaction will be formulated, Sec. [I A 



this model will used to calcu lated the dephasing rate due 
to identical TLSs, Sec. |l B| , and, finally, an average pro- 
cedure over different TLSs will be considered, Sec. H C . 



Estimates and discussion will be given in Sec. [II, while 
the conclusions will be given in Sec. |v|. 



A. Qualitative considerations 

Let us start with a toy model which illustrates the 
essence of the physics involved. Then in Sec. || the results 
will be confirmed by calculation. 

Consider the electron motion in a slowly varying po- 
tential field U(r, t). Let us calculate the phase difference 
Aip between the electron waves moving from the same 
point C along the same closed path V clockwise and 
counterclockwise, see Fig. [I]. We begin with evaluation 
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FIG. 1: A closed-loop trajectory. 



of the variation AS of the electron's action S due to the 
time variation of potential U. We assume that an elec- 
tron during its motion along the trajectory V experiences 
many scattering events against both static and dynamic 
defects, so that the trajectory can be approximated by a 
smooth curve. We have 



S 



pnds 



(1) 



where n = p/p is a unit vector parallel to the tangent to 
the curve P, ds is the length element of the curve. This 
can also be written 



ds^/2m{£ - U) . 



£ = p 2 /2m being the electron kinetic energy while m is 
the electron effective mass. Expanding this equation in 
powers of the potential energy U assumed small, one gets 

AS = - ' ,,S 



-U(s,t) = - / dtU(s t ,t) 



Here s t is the electron's coordinate on the trajectory pa- 
rameterized by time t. So U depends on time both via 
the space coordinate St and explicitly. 

Let now to be the total time of the motion of an elec- 
tron along the loop V . Accordingly, the phase variation 
in the course of a clockwise motion is 



dtU(s t ,t), 



while for the counterclockwise motion one has 



(A^)_ 



dtU(s to -t,t) . 



(2) 



(3) 



The dephasing means a non-vanishing phase difference 
= (A(p) + - (Aip)_. Thus, 



(A^=]r[(A^) 2 ± -(A^) ± (A^) 



Usin 



g Eqs. (^) and (^) one can express the above expres- 
sion through J Q to dt J Q ° dt' U(st- ,t)U(st.> ,t') where i,k — 



t, t- 



tn 



■U(s ti ,t)U{s ti ,t>) 
t. We assume that there is 



no spatial correlation between the scattering centers, 
U(st,t)U(st> ,t') tx S(st — Sf), that implies 



U(s t± ,t)U(s t >,t') oc U 2 (s,t)S(t-t') 



U(s t± ,t)U(s t ',t') tx U(s,t)U(s,t -t)S(t + t' -t ) 



Using these expressions and introducing the time corre- 
lation function of the time-dependent random potential 



U(s,t)U(s,t') = U 2 f(t-t'), U 2 = U 2 (s,t), /(0) = 1 



one obtains 



to 



(A<^) 2 cx[/ 2 / dt[l- f(2t-t )] 



If there are several mechanisms responsible for dephasing 
characterized by different coupling strengths and differ- 
ent correlation functions the resulting phase variance can 
be expressed as 



(A<p) 2 cx 



E / ~[l-/.(2*-*o)] 

„ JO T s 



(4) 



Here we have absorbed the random scattering potential 
into the partial relaxation rates T~ . They, as well as 
the correlation functions, depend on the properties of 
dynamic defects. We wish to emphasize that Eq. (Q) 
demonstrates the following point indicated above. If the 
defect has not relaxed during the time 2t — to between 
two acts of scattering then in spite of non-invariance of 
the Hamiltonian respective to the time reversal there is 
no phase relaxation. 

We distinguish two mechanisms of dephasing. The first 
is connected to interactions which cause real transitions 
between different states of the environment. This can be 
illustrated by the famous double slit experiment. If we 
send electrons at the double slit, it will pass through both 
slits and interfere with itself, creating an interference pat- 
tern on the screen. Putting detectors to determine which 
slit the electron really passed through will destroy the 
interference pattern. If the interaction with the environ- 
ment in any way allows us to determine the path of the 
electron, interference is lost. The second mechanism of 
dephasing is related to a change in the state of the en- 
vironment due to its own internal dynamics. A dynamic 
environment leads to a difference in a scattering potential 
"felt" by an electron state during clockwise and counter- 
clockwise motion. As a result, time-reversal symmetry is 
broken and the interference pattern decays. 

At this point we would like to compare our description 
to the one given in Ref. |l2] where it is proved that the 
dephasing can be described in two equivalent ways. Ei- 
ther you consider the change in the electron phase of you 
consider the change of state of the environment, where 
complete dephasing corresponds to the environment be- 
ing in orthogonal states. The last point of view would 
imply the existence of only the first mechanism of de- 
phasing that we consider. We want to emphasize that 
our second mechanism is not in conflict with this, but 
is a result of our description of the process. In Ref. [T^ 
the environment is considered as a mechanical system 
evolving according to its own Hamiltonian, whereas we 
consider the environment to be a statistical system at 
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some temperature. That is, we calculate the action of 
the environment on the electrons, but do not consider 
the action of the electrons on the environment. In prin- 
ciple, if one were to follow all the complex dynamics of 
the environment one would find that it does indeed evolve 
into orthogonal states as the electron dephases according 
to the third mechanism, and it would be seen that this 
is only the first mechanism in disguise. However, as the 
environment consists of a macroscopic number of degrees 
of freedom is it more natural to treat it statistically as a 
thermal bath. In other words, the phase of an electron 
state forming a Cooperon trajectory decays due to real 
transitions in the thermal bath formed by other electrons 
assisted by virtual processes involving dynamic defects. 
In a perturbative approach these processes occur in the 
fourth order in the electron-defect coupling constant. In 
particular, they do not enter the second-order calculation 
of defect-enhanced electron-electron interaction]^ How- 
ever, it will be shown that they play an important role 
in dephasing. 

A different classification can be made by discriminat- 
ing between the two different regimes of phase dynamics 
- phase jumps and phase wandering (or phase diffusion). 
To understand this, let us consider the first mechanism. 
Consider the case of interest for weak localization, that 
of one electron traveling around a closed loop both in the 
clockwise and counterclockwise direction, and interferes 
with itself after completing a full circuit. If we can de- 
termine which direction the electron went, we will not 
get interference. If we are unable to do so, it will ap- 
pear. As a detector we use a two level system that is 
placed at a point on the right hand side of the loop, the 
distance from the starting (and ending) point being a 
fraction a < | of the total circumference. The energy 




t=dT 



went, and the interference pattern will be smeared out, 
but not lost entirely. In this case we need the interac- 
tion with a number of two level systems along the path, 
and the combined result of all the detection times can be 
put together to determine the direction. In this case we 
speak about phase wandering or phase diffusion, since the 
random contributions of the difference two level systems 
makes the electron phase change in a diffusive way. 

Let us estimate the dephasing rates in the different 
cases. Consider first the first mechanism where we have 
inelastic electron scattering due to direct transitions be- 
tween the two TLS's states. We emphasize that the 
inelasticity is not essential to this channel of dephas- 
ing. The important point is that there is a real tran- 
sition between orthogonal states of the TLS. Since we 
are not considering degenerate states this will mean in- 
elastic scattering in our case. If the energy transfer E 
is large enough, the phase relaxation time t v is equal to 
the typical inelastic relaxation time Ti, which is a func- 
tion of the defect parameters. The criterion of "large" 
E in this case is given as Eti/H >> 1. For smaller E 
one deals with a phase diffusion or wandering. To esti- 
mate the dephasing time for this case let us recall that 
the phase coherence for any two-level system is conserved 
during the time t < h/E. While traversing the trajectory 
during the time t an electron appears to be coupled with 
N r~j t/n dynamic defects. The evolution of the electron 
wave function due to coupling with any of these defects 
is described by a phase factor exp(±iEt/K) where ± cor- 
responds to the sign of the energy transfer. If T > E the 
probabilities of the both defect states are almost equal, 
and the correlation function of the time-dependent ran- 
dom potential is f(t) — cos(Et/H), see the calculation 
later. The resulting electron phase shift turns out to 
be N 1 / 2 Et/h. Consequently, the phase relaxation time 
can be estimated as t v ~ H 2 ^ 3 t^ 3 / E 2 ^ 3 . A similar ex- 
pression for the dephasing time has been introduced in 
Refs. [l4| in connection with decoherence due to quasielas- 
tic electron-electron scattering and in Refs. [l5 16 in con- 
nection with decoherence by low- frequency phonons. In 
the following we will call this regime the phase wander- 
ing. Summarizing, we can express the contribution of 
inelastic processes as 



= max 



{^{K/E) 2 ' 3 } 



(5) 



FIG. 2: Closed loop with TLS detector. 

splitting of the two level system is E, and it starts out 
in the lower state. When the electron passes, it excites 
the two level system. We determine the direction of the 
electron by measuring the time at which this happens. 
The accuracy with which we can make this measurement 
is limited by the uncertainty principle AtAE > h. HE 
is large there is no problem, and the interference is de- 
stroyed by a single detection event. This case we call a 
phase jump. In the opposite case where E is small we can 
not determine with certainty which direction the electron 



Moving to the second mechanism, we find that the sim- 
plest way to evaluate this contribution is note that r s has 
a sense of the time at which Aip(t) « 1 provided all the 
involved defects would suffer a transition. It is clear that 
if the phase shift Sip due to transition of a single defect 
is > 1 then a single TLS is enough to produce the de- 
phasing. For Sip << 1 the significant phase evolution is 
possible only with the help of many defects. The actual 
dephasing time, t v , is also sensitive to the defect transi- 
tion rate 7. Indeed, the correlation function f(t) for sta- 
tistically independent defects is expected to have a form 



/(*) 



-2-y\t\ 



(this form will be supported by the calcu- 
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lations in Sec. |IIB| ). If 7T3 > 1, then with a help of Eq. || 
one obtains (A^) 2 = £o/ r 3 (for the reasons which will be 
clear later we ascribe the subscript 3 for the relaxation 
mechanism) . If 7T3 << 1, one has (A(p) 2 ~ 7i 2 /T 3 . We 
observe that there is a phase wandering regime also for 
this relaxation mechanism. Again, defining t v as a time 
at which (A(p) 2 ~ 1 one has 

rf=max{T3,(r 3 /7) 1/2 } • (6) 



One notes that in course of the above considerations we 
exploit the additions to the electron phase acquired by an 
electron in course of traversing of the potential induced 
by the TLSs. For each TLS the corresponding contribu- 
tion can be estimated as USr/vp, where U is the potential 
magnitude, Sr is the potential spatial scale, while vf is 
the Fermi velocity. An important note should be made 
in this concern. Since the trajectory in Eq. (jjj) and in the 
following ones is considered to be given and the positions 
of the scatterers are expected to be along the electron tra- 
jectory, the phase addition mentioned above is, strictly 
speaking, beyond the Born approximation for the elec- 
tron scattering. Indeed, in the Born approximation the 
scattering amplitude is real at least for symmetric scat- 
tering potentials. However it is possible to describe the 
phase relaxation even within the framework of the Born 
approximation if one takes into account that the "centers 
of gravity" of the two TLS states are spatially separated 
by some vector a (which is an inherent feature of the 
model suggested in Ref. [l?]) . In this case the phase vari- 
ation due to a transition within the i-th defect is simply 
given as Sip ~ (p-a)/fi. Correspondingly, the estimate for 
the proper rates in Eq. (0) is r{~3 ~ t~^(pfO,/H) 2 where 
T ed * s a typical elastic relaxation rate due to dynamic 
defectsBH 

The previous estimates are relevant to a set of defects 
having identical parameters. However, in real systems 
the defect parameters are scattered, and one has to per- 
form a proper average. As we will demonstrate, for a 
realistic model the phase wandering regime turns out to 
be important. To our knowledge, this fact has not been 
appreciated in the previous papers dealing with defect- 
induced decoherence. 

In the following sections we will give a more formal 
derivation of the dephasing rate using the Green func- 
tion method. It will permit us to consider not only the 
limiting cases but any relations between various times of 
relaxation. In the Appendix [b| we will map the results for 
the relaxation-dynamics contribution to a simple model 
of short-range defects hopping between two states sepa- 
rated some distance in real space. This model is often 
used to interpret results on the so-called random tele- 
graph noise observed in nanostructures. 



II. THEORY 

A. The model 

The dephasing mechanism is based on the assump- 
tion that in any crystalline metal there exist dynamic 
defects of a special type. These defects are tunneling 
states which are described by the Hamiltonian 

H d = {Aa 3 -Aa 1 )/2 (7) 

where A is the diagonal level splitting, A is the tunnel- 
ing amplitude, while a% are the Pauli matrices. The tun- 
neling amplitude A describes the tunneling between the 
interstitial positions while the spread of A is determined 
by (mesoscopic) disorder around the mobile defect. Con- 
sequently, we will assume that the distribution of A is 
narrow and it is centered around some value Ao. As we 
will see, one can expect smooth temperature dependence 
of dephasing at T > Ao. The above model has been pro- 
posed and successfully exploited in Ref. [l7] to interpret 
zero-bias anomalies observed in metallic point contacts. 
Note that it differs., from the well-known TLS model in 
amorphous metalstLil where the distribution of log A is as- 
sumed to be uniform, however resembles the TLS model 
for crystalline materials suggested by PhillipsE3 to de- 
scribe acoustic experiments in crystalline Si. 

Consider now spinless electrons which scatter against 
tunneling defects with the Hamiltonian (^). The total 
Hamiltonian then can be expressed in the form 

H = H d +Y, £ pCp c P + H in t (8) 
p 

where 

= \ E (^p + p' + c p + V e«P-PW» . 

pp'n 

(9) 

Here 1 is the unit matrix, V ± represent the components 
of a short-range defect potential, while r„ is the coordi- 
nate of the n-th defect. The Hamiltonian (|9j) is equivalent 
to the assumption that the electron scattering amplitudes 
are V + ± V~ in the "left" and the "right" defect posi- 
tions, respectively. Estimates for V + and V~ are given 
in Refs. |ll|,|l8[ After the transform which makes Hd di- 
agonal we arrive at the Hamiltonian 

« = ^E^3+£^p + ^E{^p + p< 

n p pp'n 

+ (|« + t«)Vp-}4ve*-'>-'f«)) 

where E n = y/ A 2 + A 2 . One observes that there are 
two processes of electron-defect interaction described by 
the items proportional to o\ and 03, respectively. They 
correspond to the two mechanisms discussed above and 
described by Eqs. (|5|) and (^). Now we proceed to more 
formal calculations in which the relaxation time t\ and 
T3 will be specified. 
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B. Quantum contribution to conductance 

The object which we will consider is the weak lo- 
calization correction to the conductivity, 5a, which for 
the case of a short-range scattering potential can be ex- 
pressed through the electron Green's functions Gr/a in 
the forrriia 



5a 



m 2 d 



(dp) (dq)p 2 



dn 
de 



ds 
2^ 



du> 
2^ 



xG R (e,p)G A (e,p) F(e,uj,p,q- p) 

xG fl (e + w,q-p)GA(£ + w,q-p). (11) 

Here d is the dimension of the problem, (dp) = 
d d p/ (2irH) d , n(e) is the Fermi function, while 
F(e,u, p, pi) is a two-particle Green's function spe- 
cific for the problem under consideration. The above 
expression describes the main contribution to the 
conductivity which arises from the region lot, qi -C 1 
where r and I are the total relaxation time and length, 
respectively. The function F(e, w, q,p) can be repre- 
sented as a sum of the maximally-crossed diagrams (the 
so-called Cooperon) which is a sum of a ladder in the 
particle-particle channel. The Cooperon satisfies the 
Dyson equation shown in Fig. [| In this figure, the 
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FIG. 3: The equation for the Cooperon. 



Cooperon is drawn as a filled square, thick lines with 
arrows correspond to the Green's functions averaged 
over the defect position, as well as over the states of the 
thermal bath, while dotted lines represent propagators 
for electron scattering against dynamic defects. Since 
the interaction Hamiltonian (|10|) contains items of 
three types, the propagator consists of a sum of three 
terms. Each propagator can be expressed as a loop 
graph where dotted lines represent Green's functions 
for a dynamic defect. To express the propagators in an 
analytical form, we will employ the technique developed 
by AbrikosovHJ According to this technique, a two-level 
system describing the dynamic defect is interpreted as a 
pseudo-Fermion particle with the Green's function 



9±(e) = (eTE/2-X + i5)- 



(12) 



where A is an auxiliary "chemical potential" which after- 
words will be tended to infinity. This trick allows one 
to remove extra unphysical states which appear since the 
Fermi operators have more extended phase space that the 
spin ones.E^Eil As a result, the Matsubara technique can 



FIG. 4: Schematic representation of the defect propagator. 



be effectively used, and after a proper analytical contin- 
uation and the limiting transition A — ► oo the quantity 
A drops out of all the expressions. As a result, the re- 
tarded propagator describing the inter-level transitions 
in the defect can be expressed as 



, E 
tanh — ; 



1 



1 



2T \w-E + iS 



E + i5 



(13) 



Here 5 is the adiabatic parameter, 5 — > +0. The propaga- 
tor describing electron-assisted transitions has the form 



i 



217 



Tcosh 2 (E/2T) tJ + 2i-f 



(14) 



Here 



j(A,E) = (j?) 70 (E), l0 (E)= X Ecoth^, (15) 

where x = 0.01 — 0.3 is dimensionless constant dependent 
on the matrix element where r yo(E) has the meaning 
of maximum hopping rate for the systems with a given 
interlevel spacing.Eil 

For the elastic component oc J we shall use a trick 
which will allow us to consider the elastic channel in 
a unified way with the inelastic ones. Namely, to keep 
proper analytical properties of the retarded Green's func- 
tion we define the elastic propagator as 



i 



i 



2T \lu + v + id lu — v + i5 



(16) 



At the final stage, the limiting case 5^0, v — > should 
be calculated. Note that the factor T _1 will be can- 
celed by the Planck function N (uj) which will appear 
in course of derivation of the equation shown in Fig. |^. 
The physical reason of this cancellation is that the elas- 
tic impurity scattering is temperature-independent. Note 
that the propagators do not include the electron-defect 
coupling constant, hence each propagator should be mul- 
tiplied by |FF«| 2 where = V+, = (A/2E)V~, 
W {3) = (A/2E)V~ . If there are additional static short- 
range defects their contribution modifies 0th propagator 
by the replacement \W^\ 2 -> \W S \ 2 + \W^\ 2 where W s 
is the contribution of static defects. 
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The equation shown in Fig. || has been analyzed follow- 
ing the procedure-af analytical continuation of Matsubara 
Green's functiont^l with making use of -analytical prop- 
erties of two-particle Green's functions £3 The resulting 
equation for F(e, ui, p, q — pi) has the form 

-F(e,^,p,q-p) = T>(u) 

M^F (£ , W ', P ', q - P )P(.-.0 
2m 

xG R (e +cj-cj', P ') G A (e + w', q - p') 
x [N (u/) - N (l>' - «)] . (17) 

Here (dp) = 2d 2 p/(2nh) 2 = pde p de/2i:, 9 is the angle 
with the x axis and we write all formulas for the most 
interesting case of a two dimensional system. 

V(u) = \W s \ 2 n s [V*(tu) - V£(wj\ 

+ ^\WW\*n d [V*{ U )-V*{ U )] . ( 18 ) 



Equation (\1T\) describes the dominant contribution pro- 
vided the sum of the incoming momenta, q, is small: 



< 1. 



(19) 



Here £ — vpr is the electron mean free path, while r is 
the electron life time, 

T- 1 = T- l +rr 1 +Tf 1 . (20) 

Here we introduce the elastic relaxation rate as a sum of 
the contributions of static and dynamic defects, t" 1 = 

T e,l + T~l With 

T~l = 2npn s \V s \ 2 /h, t~1 = 2irpn d \V+\ 2 )/h, 

and a typical inelastic relaxation rate as 

7-r 1 = 2npn d \Vj\ 2 /h, T eA /n « (p F a/K) 2 < 1 . (21) 

Here n s is the concentration of static defects while n d 
is the concentration of dynamic ones, p is the electron 
density of states. In principle we now have two sets of 
relaxation rates. From the interaction vertices we get 

r^=rr\K/2E) 2 , r^ 1 = r^(A/2E) 2 , 



while the rates appearing in (|20| ) arises in the evaluation 
of the self energy diagrams as shown in appendix A, and 
are given by 

rf 1 = rr^K^Efg^e) , = r^(A/2E) 2 g 3 (e) . 

The functions Qi^ are discussed in appendix A. We show 
that they decay exponentially at E ^> T, and only the 
regions with e < T arc important. Since we only are 
interested in this region, we will neglect the energy de- 
pendence of these functions, and put them to 1 in the 
following. The two sets of relaxation rates will then be 
the same. 



To analyze Eq. (|17]) it appears convenient to transform 
it to the form similar to the Boltzmann equation for an 
electron diffusion. For this let us take into account that 
at small q and u> the product of the Green functions in the 
integrand is a sharp function centered at e = e p i = ep 
where ep is the Fermi level. Thus it is natural to assume 
that F(e,u>, p', q — P i) depends only on q, to and the 
product q-v'. Having that in mind we first integrate over 
e p > and make use of the inequalities p F £/h ^> 1, fiw<cr 
which we assume to be met. 

The result can be expressed in terms of a new function 



F(s,uj,p,q-p) 
lu(1 — irqv) 



(22) 



where v is the electron velocity. Here we assume that 
e < T and omit the variable e. 

Following the procedure described in Ref. [l5] we ex- 
press the equation for J- in the form: 



(l + Dq 2 T)F(e,q,Lj) = 



Airpto 



-T 



diu' 



(2ni){uj - 2uj' + i/2t) 



Here D = vp£/d is the diffusion constant. Transforming 
Eq. ( p3| ) to the time representation with respect to to we 
obtain 



(1 + Dq 2 T)T(e,q,t) 



g(M) 

2tT P 



dt' 



(f-tyr 



^(e,q,t')$(£,2t-<'),(24) 



Here we denote 



^ £ ^) = - + -co S ^ + -e' 2 ^ h . (25) 
r e n n r 3 

Here the limiting transition v — > has been already 
done. Note that the function i) depends on the en- 
ergy variable e through the relaxation times r e , n and 
T3. In the following we omit the variable e in all the 
functions keeping in mind that the relaxation rates are 
energy-dependent, see Appendix A. Also in writing the 
expression for $(e,t) we have assumed that £<T. For 
E > T it decays to 1 which means that defects with 
E > T do not contribute to the dephasing, see the dis- 
cussion below Eq. ([33]) 

Equation (^4|) can be solved exactly. The solution is 
based on the relation between the kernel of the integral 
equation, 

JC{t, t') = (Dq 2 r - l)e {t '- t)/T [1 - \{2t - if)] ; 



X(t) = — ( 1 - cos ^ ) + — ( 1 - e 
n\ n I T 3 



-2-yt/h 



and its resolvent, 1Z, defined by the integral equation 
rt dt' 

—K.(tx,t')Tl(t',t) =JC(t 1 ,t)+K(t 1 ,t). (26) 
ti r 
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The relationship has the formed 

r< dt' 



^(t,q) = $(t) 



K(t,t')$(t'). (27) 



If one can construct a differential operator of the form 
Cti = J2j o,j{ti){d^ /dt{) such that 



(28) 



for any t, then the integral equation ( pTj ) is reduced to 
the differential equation (pq) for a fixed t. That can be 
directly checked applying the operator Ct 1 to relation 
(^6|). The boundary conditions corresponding for t\—*t 
can be extracted from the integral relation (Eq) and its 
derivatives with respect to t\ at t\ — > t. 

The results have the simplest form at r e <C t\, T3, 
Then one choose 



£ tl = r— + 1 



From (|2^), the differential equation for the resolvent 
1Z(t\ , t) acquires the form 

[r 3 (d 3 /dt 3 ) + T 2 {2 + X 1 )(d 2 /dt 2 ) 

+ r(l + 2Ai)(d/dt)+Ai]ft(i,0) = G(29) 

Here Ai = A + Dq 2 r. Since the phase relaxation is a slow 



process with respect to the scale r the equation ([29|) has 
small coefficients at senior derivatives which makes use- 
ful the WKB approximation. Consequently the physical 
solution can be sought in the form lZ(t, 9) cx exp[<^(t)/r] 
where <p(t) satisfies the equation, (ip + l) 2 [<p + X\(t)] — 
0. Since X(t) 1, the WKB solution corresponds to the 
equation 



<p + Xi(t) = 0. 



(30) 



The boundary condition to Eq. ( |30|) can be extracted 
from the relation lZ(t,t) = 0. In this way, we obtain the 
quasiclassical solution in the form 



1Z(t, ti) = exp 



* dt' 



Ai(t': 



(31) 



Now we substitute Eq. (|3f|) in the expression (J27j) to ob- 
tain the final expression for the Cooperon T . The first 
item in Eq. (|2^) is the contribution of lowest order scat- 
tering and it should be neglected in the diffusion approx- 
imation. Here we analyze the quantum contribution to 
the static conductance, so only .F(O) is important. As a 
result, we obtain 



^(0,q) 

m 



f° ^(tVCO, (32) 

J-oo T 

sm(Et/h)~ 

En/h 

(e 2 ^ t/h -l) , << 0.(33) 



Dq 2 t 



t 

n 



t 

T3 



27T3 



An important feature of Eq. (|33j) is that if one neglects 
the processes in which the defect changes its state then 
the dephasing is absent. Indeed, putting t\ = T3 = 00 we 
get = 1 and T{q) = (Dq 2 )- 1 . This results in log- 
arithmic divergence of the conductance in the 2D case. 
Another important feature is that at small time t, which 
has a physical meaning of the time difference for the 
collision act for clockwise and counter-clockwise partial 
waves, no linear in t term is originated by inelastic pro- 
cesses. Physically it means that no dephasing takes place 
if the scattering defect had no enough time to change its 
state. Therefore, the dephasing appears proportional to 
the probability for the defect to escape the state in which 
it has been registered by one partial wave. 

Substitution Eq. (||) into Eq. (||) and then into 
Eq. ([ll]) we obtain for 2D case the expression for the 
quantum contribution to the conductance. Since only 
l e l T -C £f arc important one can neglect e- 
dependence of the relaxation rates and put e = in the 
expressions for the these quantities. As a result, one ar- 
rives at the well-known expression, 



2Tr 2 h 



ln^ 



where t v is defined according to the equation 



ln^ 

T 



'dr] 
— i 

1 V 



Ti 



sin(?7 Et/K) 



(34) 



Et/K 



T 3 (ri,E,A) = - 

T3 



1 2 7 t v ' 



where rj = i/r. This equation is obtained by the integra- 
tion over q. 

We can now recover our estimates (^) and (|^) by the 
approximate value of the integral 

1 V 

in the case where a <C 1. We split the integral at the 
point rj*: 



I = 



n 



v 



Defining 77* such that mf n — 1, we have 77* > 1. The 
integrand will be very small in the last integral, and we 
neglect this. In the first integral we put the exponent 
equal to 0, and get 



I « In rj* = In a 



-l/r. 



(35) 



Expression (|3J) depends upon two dimensionless quan- 
tities. The first one is Et/h. As follows from Eq. (pq 1 ), 
the typical value of t is t v . This parameter determines 
the efficiency of the first mechanism of dephasing, that of 
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direct transitions of the TLS states. If Er^/h < lwe are 
in the regime of phase wandering, and we can expand I\ 
in Eq. (Q) in powers of this parameter up to the lowest 
order. However, if Et^/H > 1 we have the case of phase 
jumps, and we can neglect the sine term in Ti. These 
expansions are given in the equations ([58|) and (^) be- 
low. In both cases we can use the formula (^||) to arrive 
at the estimate (^). 

The second dimensionless parameter is jt^/H. It de- 
scribes the effect of the second mechanism of dephas- 
ing arising from the 173 vertex. The physical expla- 
nation is that the dephasing occurs only if the partial 
waves meet the scatterer in different states. Expand- 
ing in small [phase wandering, Eq. p^)] and large [phase 
jumps, Eq. (46J)] values of this parameter we arrive at the 
estimates (|). 

In addition there is also the dimensionless parameter 
7T1//1 which will control the effect of the second mecha- 
nism acting through the o~± vertex. This has been ne- 
glected in the above calculations since the inequality 
7 <C E/h is met. 

If the estimates (Q) and (|6|) have different orders of 
magnitude then the shortest one is effective. However, 

l/r v + 1/t« + 1/rf 

since Ti and 1?3 depend on time in different ways. The 
most clear manifestation of this fact is seen in the mag- 
netic field dependence of the quantum contribution. 



C. Average over different dynamic defects 

To calculate the quantum contribution to conductance 
one has to sum over different dynamic defects. In the pre- 
vious considerations we have assumed that all dynamic 
defect have the same interlevel distance, E, and the same 
hopping rate, 7. Consequently, the summation over dif- 
ferent defects has been allowed for by the factor na in 
the expressions for the relaxation times Tj. However, in 
realistic systems both E and 7 can be distributed over a 
significant range. Since the number of dynamic defects 
at a typical electron trajectory is assumed to be large 
the summation over different defects can be replaced by 
a proper average. To calculate the latter it is necessary 
to specify the distribution function r P{E, r y) which we as- 
sume to be normalized to 1. To specify this function, 
let us come back to the effective Hamiltonian (^). Since 
A is determined by the defect's neighborhood while A 
is determined by the distance between two metastable 
states it is natural to assume A and A to be uncorre- 
cted, V(A,A) = V a (A)Va{A). 

Below we will discriminate between two model distri- 
butions. The fixftUone will be referred to as the "glass- 
model" (CM) Hid According to this model the distri- 
bution of A is assumed to be smooth, V\ = Vo- Since 
the tunneling integral A is an exponential function of 
the distance between the potential minima and the lat- 
ter is smoothly distributed, it is assumed that V\ oc A . 



Within this model it is natural to choose the interlevel 
splitting E and the quantity p = (A/E) 2 as independent 
parameters. Since 7 oc p, equation ( J15| ) can be rewrit- 
ten as 7 = pjo(E). Consequently, the GM results in the 
exponentially-broad distribution of relaxation rates. Fur- 
thermore, to keep the distribution normalized, we intro- 
duce cut-off Pmin(-E) = 7min (E)/jo (E) <C 1 and assume 
C = m(l/p min ) = ln(7o/7 m in) to be energy-independent. 
A cut-off energy in the smooth distribution of E at some 
E* is also assumed. As a result, we get the distribution 



Vgm(E, P ) 



Q(E* - E) 1 



E*C p^T=p 



(36) 



Another model which we will call the "tunncling- 
states-model" (TM) is more appropriate for crystalline 
materials. There the tunneling integrals A is determined 
by the crystalline structure and are almost the same for 
all dynamical defects. On the other hand, the parameter 
A is determined by long-range interactions, and it is as- 
sumed distributed smoothly within some band, cf. with 
Rcf. [L7|. Then 

V TM (E, A) = @iE *~ E) -f=£= 6(A - A ) . (37) 



E* 



A 2 



In the following we will assume that the dynamical de- 
fects are characterized by Ao < T. To calculate the 
total contribution of the dynamical defects in the case 
when their parameters are random one has to replace I\ 
in Eq. @) by the averages f 2 (ry) = / dE dA r,(ry, E, A). 
Below we will discuss in detail only the tunneling-states- 
model which seems to be more appropriate for crystalline 
materials. 

Let us discuss the contribution of direct transitions and 
relaxation separately. 

a. Contribution of direct (resonant) transitions. 
The item Ti, responsible for the direct transitions, 
is proportional to |J-Fi| 2 = (A/£) 2 |U (1) | 2 . Since 
Vtm(E,A) oc (5 (A — Ao) the integral over A yields the 

factor A 2 /E^/E 2 - A 2 . Using @, the quantity t£ 1} is 
estimated from the expression f i (i~!p /rj = 1. The fol- 
lowing calculation depends on the relationship between 
E and t v . At Et^/Yi <C 1 one can expand the expression 
for Fi as 



T 1 (r ] ,E,A)^(A/h) 2 (r,T) 3 /T l 



while at Er v /h ^> 1 



T 1 (r ) ,E,A)~(A/E) 2 r ) T/T i 



(38) 



(39) 



To estimate Yi{rj) let us introduce the energy splitting 
Ta at which Er v /h = 1. The meaning of this is that 
a TLS with an energy splitting less than T\ will not by 
itself cause complete phase loss, i.e. we are in the regime 
of phase wandering. A TLS with E > T\ causes a phase 
jump. Let us first assume that 



A <Ta<T. 



(40) 
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Using the distribution (|37]) and returning to dimensional 
time one obtains 



Fi(t) 



t 



A 2 
T A E* 



t 



h 



Ta 
E* 



(41) 



Now, let us define r^p and Ta to make both contributions 
to r 1 (r^ 1 ^) equal to 1. This definition of T A is consistent 
with that given above within the accuracy of the approxi- 
mation because the two terms are the expansions in large 
and small values of Er^/h. The point where both terms 
becomes of the order 1 should then correspond to the 
crossover point Et^/H = 1. This is easily checked from 
the formulas below. In this way we get 



rW=r A (T A /Ao) 



where 



T A = (Mq/ta) 1 / 2 r A = n(E*/A ) . 



(42) 



(43) 



The time ta is due to the dynamic defects with symmetric 
potentials and energy splitting equal to A . At T < 
T A for all energies the inequality Et^/H -C 1 is met, 
and only the second item in Eq, pl| ) is important. One 

should replace Ta by T in this expression to obtain — 
t a (T a /A )(T a /T) 1 / 3 . In contrast, if T A < A then only 

the first item in Eq, ([ll]) is important. In this case rip = 

TA- 

The result can be summarized as 



1 



l_Ao 
r A T A 1 T A /A 



^(T/Ta) 1 / 3 ,!}, T a » A , 



T a « A c 



(44) 



b. Contribution of relaxation processes. Since only 
E < T are important, for estimates one can assume 
EcothE/2T « 2T. Thus, 7o w 2 X T becomes E- 
independent. For the same reason t^ 1 can be approx- 
imated as t~ 1 (A/E). The following calculation depends 
on the relationship between 7 and t v . At jt v -C 1 one 
can expand the expression for 1^ as 

T 3 ( V ,E,A)^r] 2 7 T 2 /r 3 = 7o ( V r) 2 T^A 2 A 2 /E i . (45) 

while at jt v ^> 1 one has 

T 3 = V(r/r 3 ) cx (A/E) 2 . (46) 

To estimate 1^3(77) let us introduce the energy splitting 
E x at which jt v = 1. The meaning of this is that a TLS 
with E < E x will probably jump during the trajectory 
traversal time t v (it is "fast moving"), whereas a TLS 
with E > E x will have a low probability to jump in the 
same time (it is "slow moving"). First we assume that 



A <£v<CT. 



(47) 



Using the distribution ( |37| ) and returning to dimensional 
time one obtains 



TV E* 



hn E X E* 



(48) 



Now, let us define t v and E x to make both contribu- 
tions to f 3 equal to 1. This procedure indicates that the 
defects with E = E x are those which experience a hop 
during the typical Cooperon trajectory traversal time. 
One obtains 



V4 3) = rr\E x /E*) , E x = A { X Tr A /h) 



1/3 



(49) 



Introducing the characteristic temperature T a — H/xt a 
at which E x — Aq one can express the dephasing rate as 



(3) _ 



-\T/T a ) 



1/3 



(50) 



Another important characteristic energy is the temper- 
ature Tp at which E x = T, Tp = A 3/2 /Ti /2 The ratio 
Ta/Tp = (xA t a /H) 3 / 2 can be arbitrary. 

The meaning of T a and Tp can be understood as fol- 
lows. Imagine starting at some large temperature where 
T 3> E x ^> Ao. As we lower the temperature E x is 
also decreasing, but it decreases at a slower rate than 
T (E x - T 1 / 3 ). At the temperature T Q , T; x = A . 
Since Ao is the lower cutoff for E, if T < T Q all defects 
are slow-moving because no defects exist with sufficiently 
low splitting. Alternatively, since T is decreasing faster 
than E x , T will overtake E x at the temperature Tp. For 
T < Tp all defects are fast, because the slow ones are 
frozen out. Which temperature is reached first depends 
on the specific values of the parameters (since the ratio 
T a /Tp is arbitrary). Also, if T a > Tp then E x < A 
for all T < T a , so in particular Tp < Aq and is thus 
unimportant. Similarly if T a < Tp then T a < Aq. 

The result @ is valid if T > E x > A , or at T > 
T ai Tp. At E x > T > A , or Tp > T > A , only the 
first item in Eq. ( E ) exists and the quantity S x should 
be replaced by T. As a result, 



l/r« 



-'(T/Ao). 



(51) 



Since A = Tp /3 T^ /3 the results @ and <^f match at 
T = Tp This temperature region exists only if Tp > T a . 
At T < A the dephasing rate strongly decreases with 
the temperature decrease. 

If T a ^> T ^> Tp the relaxation is slow for all the dy- 
namic defects and the second item in Eq. ( [48| ) is impor- 
tant. However, in this case one has to replace E x — > Ao 
in its estimate. As a result, 



- 1 (T/T Q ) 1/2 



(52) 



The temperature dependence of t v is sketched in Fig. || 
for T Q < Tp and in Fig. | for T a > Tp. 

c. Comparison Now we are in a position to compare 
the contributions to the dephasing rate. Both contribu- 
tions to 1/t v are parametrized by the quantity 1/ta, the 
relative contributions being dependent on the temper- 
ature. The relative resonant contribution crosses over 
from (A /T A )(T/T A ) 1 / 3 to A /T A at T = T A . The rel- 
ative relaxation contribution crosses over from T/Aq to 
(T/T a ) 1/3 at T = Tp if Tp >T a . In the opposite case 
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In T 




FIG. 5: Schematic picture of In — as function of temperature 
for T a < Tp. 



In x 



interesting regimes it is the quantity A 2 that has to be 
averaged. The results then can be expressed as 



AT" 



d£,e 



at An ^S> A 



As it is seen, the only effect of the scatter in A corresponds 
to renormalization of the tunneling matrix clement by a 
constant factor e . Since this factor is temperature- 
independent, it does not affect qualitatively the picture 
obtained with an assumption of a fixed value of A. One 
also notes that the averaging procedure practically cuts 
out a contribution of the region A > Ao since at this re- 
gion the tunneling matrix element exponentially decays 
with A. Thus the picture is not sensitive to this region, 
and any distribution of A with a lower cut-off (even allow- 
ing a Gaussian smearing of this cut-off) does not change 
the considerations made in an assumption of a fixed value 
of A = An. 



III. ESTIMATES AND DISCUSSION 



FIG. 6: Schematic picture of In as function of temperature 
for T a >T(3. 

is crosses over to (T/Tq,) 1 / 2 at T = Ao and then to 
(T/T a y/ 3 at T = T a . 

We conclude that at T > T a ,Tp the relaxation contri- 
bution dominates, and the dephasing rate is 



l(T/T 



a/3 



•C] 



(53) 



where £ is a constant of the order 1 originating from the 
resonant contribution. 

At low temperatures both contributions can be impor- 
tant, their interplay depending on the relationship be- 
tween the temperature T and the characteristic energies 
Ao, Ta, T a and Tp. For both mechanisms the dephasing 
rate vanishes as T — > and there is is a region in which 
the dephasing rate is proportional to T 1 / 3 . 

d. Averaging over the tunneling matrix element. In 
our considerations we have assumed that Vtm °c S(A — 
Ao), i. e. that the tunneling matrix element A is given. 
Note, however, that due to disorder the barrier param- 
eters are also scattered. To discuss role of such a scat- 
ter let us assume that the overlap integral A is given 
by the expression A = (Jtuq/tt) e _A where ujq is some 
attempt frequency while the barrier strength A is dis- 
tributed according to Gaussian law around some central 
value Ao = In ?kJo/7rAo 1, 



(54) 



Here M is a proper normalization factor which at A <C 
Aq is equal to (27rA 2 ) -1 / 2 . As we have seen, for most 



To make estimates we rewrite the expression (21) for 
r, in the form 



T; = a in v F n d ■ 



d 



<\v-/v^ 



(55) 



where <7g is the cross-section of elastic electron scattering 
by a dynamic defect. Correspondingly, the key parameter 
of our theory, ta , is given as 



t a 1 = A a P d a in v F 



(56) 



where Pd = rid/ E* is the density of states of the dynamic 
defects. 

The density of states Pd can be, in principle, estimated 
for a given material on the basis of point contact mea- 
surements. Namely, metallic point contacts are known 
to exhibit, first, telegraph resistance noiseE3 and, second, 
zero-bias anomalies;c3 both effects axe-expected to be as- 
sociated with the dynamic defects. OEffij 

Although we appreciate that the material prepara- 
tion procedure can significantly affect the defect system, 
we believe that such experiments can provide more or 
less reasonable estimates for Pd. The telegraph noise 
studiesc3 for a Co nanoconstriction with a size of ~ 10 nm 
revealed the presence of about several dynamic defects at 
energies less than 10 mV. This would give us the value 
Pd ~ (3 — 5) x 10 32 erg _1 cm —1 . However, the telegraph 
noise is related to TLS with rather slow relaxation rates 
(< 10 3 s _1 ) while we are interested in the defects with 
switching times of the order of 10~ 9 s. Consequently, 
these estimates most probably significantly underesti- 
mate Pd- What is more instructive, the magnitude of 
the resistance noise revealed rather large defect asymme- 
try corresponding to the estimate a- m ~ ~ 10~ 15 cm 2 . 

We believe that the zero bias anomalies can give 
more reliable information concerning Pd- The magni- 
tude of these anomalies for Co nanoconstrictionscj of 
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10 



34 



the same type as mentioned above corresponds to a 
presence of se*acal tens of TLS at the energy region 
about 1 meV.tllEa Correspondingly, one obtains Pd ~ 
(3 - 5) x 10 34 erg- 1 cm- 3 . 

Based on these estimates and taking Pd 
erg -1 cm~ 3 , <r m « 1CP 15 cm 2 , vf ~ 10 8 cm/s, and 
Ao ~ 10 mK we obtain t\ w 10~ 9 s. Equations (pq) and 
( |56| ) yield Ta — Ao. Thus at temperatures larger than 
T\ w A w 10 mK one expects, according to Eq. (44), 
temperature-independent contribution of resonant pro- 
cesses. 

For the relaxation channel, one obtains T a « Xjg w 10 
mK. Consequently, at T > T a «_7a ~ 10 mK one expects 

10~ 9 s. 



that dephasing rate obeys Eq. fl53j ) with ta 
Now let us check if our assumption Ao 
realistic. We will exploit a crude estimate 



hu (2 

A exp -- 

7r \ n 



dr v / 2MU(r) 



where U(r) is a potential relief between the two sta- 
ble defect positions separated by a distance a, and M 
is the defect mass. Taking as an example U(r) — 
(Uq/2) [1 — cos(27rr/a)] one obtains for the exponent 
(2a/Trh) y / 2U M. Taking for a light defect uj ~ 10 14 
s -1 and assuming a w 10~ 8 cm, Uq ~ 0.2 eV one es- 
timates that the value A = 10 mK is achievable for 
M w 2 x 10 -23 g which corresponds to atomic weight 
w 10. 

Summarizing our estimates, we can conclude that for 
realistic parameters of the dynamic defects one can in- 
deed expect a slow temperature dependence of the de- 
phasing rate given by Eq. (|55j ) crossing over to a rapid de- 
crease at low temperatures. The crossover temperature, 
as well as the behavior below that temperature, depends 
on the distribution of A. For a delta-like distribution ( p7j ) 
the TLS spectrum has a gap of A . Thus the TLS con- 
tribution to dephasing rate is exponentially frozen out at 
for T < Ao, and we are left with the "standard" mech- 
anisms like electron-electron scattering. However for the 
Gaussian distribution (|54|) with A ^> 1 the situation is 
different. In this case the cut-off temperature is given 
by the renormalized tunneling coupling, Aoe A while for 
lower temperatures one deals with rather flat distribu- 
tion of A within the region A < Ao + A. Correspondingly, 
at these temperatures one deals with a glass-like TLS 
distribution for which t v oc T -1 . 

Although some papers, e. g. Refs. [|^7| stated that 
to explain the dephasing saturation by a TLS contribu- 
tion one would need an unreasonably large concentra- 
tion of the TLS, this conclusion was mainly based on the 
"glassy" model of the TLS while we exploited the tun- 
neling state model of Refs. [l7],|lj| In general, to obtain 
independent information concerning the TLS concentra- 
tion based on "bulk" measurements like acoustic mea- 
surements or heat capacity measurements in conductors 
is rather difficult due to a presence of electronic contri- 
butions. In particular, the value Pd ~ 10 34 erg" 1 cm -3 
exploited above is still less than the electron density of 



states (~ 10 erg - cm -3 for Co) and so the TLS are 
not expected to affect significantly the properties of the 
material, like heat capacity. 

Now we would like to compare our results with the 
calculations given in Refs. [?]]|,^| where a similar problem 
was considered. The authors of Ref. [l8] gave a semi- 
phenomenological treatment of the problem. They ex- 
ploited the TLS distribution typical for the standard 
glassy TLS model, but with the upper cut-off Ao jmax 
for the tunneling matrix element. For the resulting de- 
phasing time they reported a proportionality of r -1 to 
Ao,max/S* (in our notations) in the limit T v Ao imax > h. 
One notes that such a proportionality is in agreement 
with the second line of our Eq. (p4h although the to- 



10 mK is tal expression for t v was some different from ours. Fur- 
thermore, the estimate for the opposite limiting case was 
completely different from our Eq. (|44[). Since Eq. ( |4^ ) 
(57) corresponds to the "resonant" or "inelastic" channel we 
conclude that the authors of Ref. [l^ accounted for only 
these inelastic processes of electron dephasing. The "elas- 
tic" , or 03 , channel (which, as we have seen, can dominate 
with respect to the "inelastic" one) seems to stay beyond 
the quantitative results of Ref. |l8|. 

In Ref. H the dephasing due to dynamical defects was 
treated within the framework of the two-channel Kondo- 
model. We believe that this model is not relevant to the 
metallic samples we are interested in, see Ref. |2§| . 

Then, the dephasing by TLS was also considered in the 
recent paper Ref. |t] where the saturation behavior of t v 
in quantum dots for the TLS distribution with fixed Ao 
was claimed. As it follows from the derivation)?] only the 
transitions between TLS states due to interaction with 
the electron forming the interference loop are taken into 
account. At the same time, the transitions due to other 
electrons forming a thermal bath (second mechanism of 
dephasing) are ignored. Hence, only the a\ channel is 
taken into account, and the result is similar to the second 
item in our Eq. (|49|). However, the main contribution 
arising from the 03 channel is omitted. 

It is worthwhile to mention that similar ideas_svere used 
to explain the magnetoresistance of polymers.Ej Polymer 
systems exhibiting rather large fraction of free volume are 
expected to form readily mobile and metastable defects 
of different types. 



IV. CONCLUSIONS 

To conclude, we have shown that the dynamic defects 
can be responsible for the slowing down of the temper- 
ature dependence of the dephasing rate at low tempera- 
tures. 

There are two mechanisms of dephasing. The first one 
corresponds to direct inelastic scattering of electrons by 
the defects, while the second one is due to violation of 
the time reversal symmetry caused by fluctuations of the 
scattering potential. The first mechanism can indeed lead 
to the saturation, while the second one still contains a 



12 



temperature dependence although a weak one. 

However, when T < A the dcphasing rate rapidly 
tends to 0. 
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Using Eq. ( |l3| ) we obtain 

Qx = n(e + E) n(e - E) rC 2 (e) . (A5) 

At E 3> T it is proportional to e~ E l T at any finite e. 

While calculating Q$ one can expand N(ui) ~ T/ui, 
n(e — oj)/n(e) w 1. After that 

G 3 = cosh~ 2 (E/2T) . (A6) 

It can be easily shown that Qq = 1. 



APPENDIX A: CALCULATION OF THE 
RELAXATION RATES 

The relaxation rate is determined as an imaginary part 
of the analytically-continued Matsubara self-energy. In 
general, we have 3 contribution to the self energy due 
to three different types of the electron-TLS interaction. 
Since for a short-range scattering potential the interac- 
tion vertexes do not have an internal structure, for each 
bosonic propagator T>i one obtains: 

Sf(e fe ) = pg 2 J d^FtHsk,^), 

F t M {e k ,i P ) = Tj^Vf(u; s )G(e k -u; s ,C P ). (Al) 



Here p is the electron density of states, lo s = 2irsT, 
£fc = 27r(fc+l)T, £ p = p 2 /2m—p, and gi is the proper cou- 
pling constant determined by the Hamiltonian (JlOj) . The 
analytical continuation is performed in a usual way. Since 
there are two cuts in the complex w-plane, at 3w = 
and 9 (e — uo) = for each i we get F R (et, £ p ) = i*i + F 2 
where 



Fx = 
F 2 = 



00 + £ fc 



N{uj)V R (uj)duj 



-00+£ fc 

oo 



2ni 

N(u))G R (u>)du> 
2wi 



[G A {e k -w)- G B .{£k - w)] 
V r (lu)-V a {u)] . (A2) 



For brevity we omit the arguments £ p pf the electron 
Green's functions. Now we replace the integration vari- 
able in the expression for Ei as w — > uj — e k and 
then combine two integrals. Taking in account that 
N (uj + iirT) = -n(u) where n(ui) = (e w / T + l)~\ 
9?Gfl(e, £ p ) — irS(e — £ p ) and making a straightforward 
algebra we obtain 



coth 



(-) 

\2TJ 



tanh 



2T 



x 3 \p R {oj)] 8{e-w- Q . 

Performing trivial integration over £ p and taking into ac- 
count the 1/2t = — $jE# we finally obtain 

Tr\e) = 2irpg 2 g i (e), (A3) 

p OO 

ft(e) = / dwJV(w)n(e-w)n -1 (e)9f [X>f (w)](A4) 



APPENDIX B: MAPPING TO A 
RANDOM-TELEGRAPH-NOISE MODEL 

Consider an electron trajectory with the total traversal 
time tg which contains N dynamic defects able to hop 
between two sites. They are rather rare, so a typical 
neighbor of any active dynamic defect is a static one . 
The total length of the trajectory is 



C = t/v f 



M 

E 

s=l 



where M is the total number of defects, M 3> N. 

Let us parameterize the electron motion along the tra- 
jectory by time t and allow some of the defects (labeled 
by j) to make transitions between their states. For those 
defects, 

Rj(t) =B,f ) +u j (t). 

The length of distorted trajectory C + traversed in the 
positive direction is 



M 



E *S 

s=l 
£(0) 



Ri 0) -u s (t s ) 



N 



I) 1 > V,- • Uj(tj 



/ , '3 
3=1 



Here = (vj_i — \j) is the change in the electron ve- 
locity due to scattering by jth EF. 

Now we can specify the displacement of jth EF as 

a,-(f) = aj£j(t) 

where £j(t) is a random telegraph process (RTP), i.e. 
a function switching between the values ±1 at random 
times and having the correlation function 



,-2 7i |t-t'| 



Then the time-dependent contribution to the length is 

(SC)j (t) = lj (t) , lj EE ( V j • Bj)/V . 

For a given defect j, the phase difference is just 
(S^) j (r) = (p F l j /h) WA-tQto-tj)]. 
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Let us split the calculation of the average (e l<5 *) in two 
steps. First let us calculate the average over different 
realizations of a given RTP, 

k(r) = (e iJ ^-^- t A , J = PF l/h. 

\ I RTP 

This sum can be calculated using the generation function 
(for tp>t a ) 

K(x,y) = (e-***^-**^ 

= e _7(t ' 3 -*° ) [cos(x + y) cosb-7(t/3 - t a ) 
+ cos(x — y) sinh y(tp — t a )] . 

Substituting x = —y = J and considering arbitrary times 
we obtain 

k(t, t ) = 2 cos 2 J + 2 sin 2 J e -^\to-2t\ _ 

We observe that the function depends explicitly on the 
position of the scatterer along the trajectory, that is nat- 
ural. It does not contain complete destruction of the in- 
terference - it appears only after averaging over different 
dynamic defects. 

The average over different dynamic defects will be per- 
formed using the Holtsmark procedure (see, e. g., Ref.^ 
for a review) according to which 

(e^) d = e- w ^ , W(t ) = n cS V c n(to) 

Here n e g is the concentration of "active" defects, V c is 
the "contact volume" , while 

W(t Q ) = (1 - k(t,t )) d = (2 sin 2 J V {t - 2t)) d 



where T](i) = 1 — e 27 '*'. The contact volume is estimated 
as V c = erupted where a is the scattering cross section. 

Let us for simplicity assume that the hopping distances 
dj of the defects are the same. We have also to assume 
that the phase changes due to individual hops to be small 
to keep the treatment consistent. Assuming J < 1 wc 
easily average over the directions of hops to get J 2 = 
(47t 2 /3)(pf«/^) 2 . Now let us average over the positions 
of the defects along the trajectories. This is done as 



K(t ) = 2J 2 / V (t - t) ^ = (2 7 - 1 + e" 2 ^) . 
Jo to ^7 



We observe that n(to) ~ J 2 min{7io, 1}. Collecting the 
factors, we obtain 



W(t ) Pa (t/r 3 ) min{ 7 t , 1} 



where 



r 3 1 « (4:TT 2 n c ff/3)(pFa/h) 2 avF . 



The concentration n e ff depends on the distribution of the 
TLS parameters. After evaluating it in a proper way one 
recovers the results of Eq. (fid]). 



* Electronic address: joakim.bergli@fys.uio.no 
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